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Stochastic noise estimator method is a powerful tool to calculate the disconnected insertion involving quark 
loops. We study the variance reduction technique with unbiased subtraction. We use the complex Z2 noise to 
calculate the quark loops on a 16 3 x 24 lattice with (3 — 6.0 and re = 0.154. Unbiased subtraction method is 
performed by using hoping parameter expansion. We report on the variance reduction for the point-split vector 
current as a function of the number of subtraction terms and the number of noise used. 



1. INTRODUCTION 

Disconnected quark loops originates from the 
vacuum polarization of sea quarks of the inter- 
acting operator. Whether we do lattice calcula- 
tion for quenched or dynamical QCD, these quark 
loops contribute significantly to various hadronic 
matrix elements. For example, for nucleonic 
scalar (irNN-a term) vector (strange mag- 
netic moment) [^|]3j, flavor-singlet axial (quark 
spin) 0, and tensor (quark orbital angular mo- 
mentum) H channels, it is absolutely essential to 
consider disconnected quark loops. In particular, 
the strangeness content of the nucleon comes ex- 
clusively from these disconnected loops. 

However, disconnected quark loops are difficult 
to simulate as they contain both diagonal and off- 
diagonal elements of the large inverse fermion ma- 
trix (M). Even for a moderate size lattice (e.g. 
16 3 x 24) one needs to invert a 10 6 x 10 6 matrix 
which requires an enormous amount of compu- 
tation time. Instead, one can use noise method 
to estimate required traces [HQ to obtain discon- 
nected loop contribution. This study focuses on 
the stochastic estimation of the strangeness mag- 
netic form factor (G S M (0)) of the nucleon and re- 
ports a systematic analysis about the number of 
noises required to extract its signal. 

Experimentally we are still uncertain about the 
value (even the sign) of the strangeness mag- 
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netic moment of the nucleon. SAMPLE jj| and 
HAPPEX expts. reported G" M (0) = 0.01 ± 
0.29 ± 0.31 ± 0.07 and G% + 0.39G^(atQ 2 = 
0.477GeU 2 ) = 0.025 ±0.020 ±0.014, respectively. 
Prediction from theoretical models vary in a wide 
range (-0.75 to +0.30fi N ) @. Lattice QCD re- 
sults also differ in conclusion. Our previous stud- 
ies § suggest G S M (0) = -0.28 ± 0.10, while ref. 
H reported very tiny signal for G S M . This work is 
aimed at studying why there is such a difference. 



2. Noise method and unbiased subtraction 

Disconnected quark loop calculation by 
stochastic noise method have been detailed in 
refs. [||,[l],|H,i|,|§ • They used random noises to 
estimate various traces involving fermion matrix. 
Using a set of random noise vectors rj, one can 
estimate the trace of a N x N matrix A as [§|,0J|] 



Tr(.4) 



E [< rfAri >] 



(1) 



Of course, this trace will be an approximation for 
finite number of noises and the variance of the 
estimator depends on the choice of the noise. It 
has been demonstrated that Z2 noise [pf is the 
optimal noise with minimum variance 0. For 
a given L number of Z2 noises, variance of this 
estimation is given by | 
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This variance can further be reduced by the 
method of unbiased subtraction [p|, where a set 
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of P traceless matrices (Q) are subtracted from 
the matrix A as 



Tt(j4) = E 



\ p=l ■> 



(3) 



where A's are some variational coefficients. Cor- 
responding reduced variance will be M 
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This subtraction is unbiased in the sense that it 
does not change the expectation value of Tr(^4). 
Choice of these traceless subtraction matrices (Q) 
should be such that they match the off-diagonal 
behavior of the matrix A. For disconnected loop 
calculation A will be replaced by the inverse 
fermion matrix M . Previously it was shown 
that the above variance for M _1 can be reduced 
substantially [|]]||||| by using a set of traceless 
matrices obtained from the hoping parameter ex- 
pansion of the fermion matrix M as 



Ar 1 = I + kD 



k 2 D 2 



k 3 D 3 + n i D i 



For the point split conserved current, 
nected quark loop can be written as 



(5) 
discon- 
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Before subtracting each matrix (I, kD, k 2 D 2 etc. 
of Eq.(5)) from M _1 in Eq.(6), one should make 
sure that it does not change the loop expectation 
value. In fact, all matrices with M _1 substituted 
with even order of D are traceless in Eq.(6). First 
and second terms (/ and kD) are also traceless. 
However, starting from k 3 D 3 , all odd orders in 
D are not traceless. So, to subtract an odd order 
term another matrix is need to be subtracted form 
it so that the resulting matrix is traceless. For 
example, for k 3 D 3 term, one needs to subtract 
following plaquette terms from the loop: 



8k 3 ^ E 
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corresponding to current in the (1 + 7^) and 
(1 —7^) respectively. One should notice that in 
each term plaquettes are at position x and x — v, 
and so, at \q\ ^ one cannot use translational 
invariance due to the Fourier transformation fac- 
tor. Similarly, one needs to subtract some chair 
diagrams from k 5 D 5 terms to make it traceless. 




x + v 



Figure 1. Plaquette term associated with current 
corresponding to 1 + 7 M . This needs to be sub- 
tracted from k 3 D 3 term to make it traceless. 

3. Results 

Numerical simulation was done on a 16 3 x 24 
lattice at k — 0.154 with 60 configurations where 
each configuration is separated by 20,000 sweeps. 
Unbiased subtraction is done with terms up to 
k 4 D 4 . We systematically study the signal for the 
strange quark form factor as a function of the 
number of complex Z 2 noises used per configura- 
tion. In Fig. 2 we plot the summed ratio of three 
to two point functions as a function of the time 
slice, from which one can obtain the magnetic 
form factor (see ref Q for notations) . Valence and 
sea quark mass is kept fixed at k = 0.154. First 
sub-figure is with 300 noises without any unbi- 
ased subtraction. Next 4 sub-figures are results 
with unbiased subtraction with different number 
of noises (30, 100, 200 and 300, respectively). It 
is clear from these sub-figures that we do not find 
any signal up to 200 noises and the signal becomes 
prominent at around 300 noises. The fitted slopes 
for 100, 200, and 300 noises are -0.052 ± 0.09, 
-0.060 ± 0.048 and -0.092 ± 0.040, respectively. 
This implies that the signal can only be extracted 
out at around 300 Z2 noises. Slope for the 300 
noise case agrees well to our previous calculation 
where we used subtraction terms up to k 2 D 2 . 
Since this result and previous result agree at one 
k, we do not carry out calculation for other n 
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values. In our previous calculation we obtained 
G" M (0) = -0.28 ±0.10 and this systemic study of 
noise versus signal supports that result. 

4. SUMMARY 

We use noise method to extract the discon- 
nected quark loops. As an example, we choose 
the strangeness magnetic form factor of the nu- 
cleon. An unbiased subtraction method is em- 
ployed to reduce the variance in trace estimation. 
This study suggests that certain minimum num- 
ber of Z 2 noises are required to extract the signal. 
In the case of the strangeness magnetic form fac- 
tor we need around 300 complex Z 2 noises. Re- 
sults of this study is consistent with our previous 
results We believe, this also explains why 
with 60 real Z 2 noises, the work of [|] did not see 
a signal even with a larger number of gauge con- 
figurations. In future we hope to carry out this 
strangeness calculation with the overlap fermion. 
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Figure 2. Summed ratio of three to two point 
functions for different number of noises from 30 
to 300. m is the fitted slope which is related to 
the magnetic form factor 0. 



